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Electron capture and /3 decay play important roles in the evolution of pre-supernovae stars and 
their eventual core collapse. These rates are normally predicted through shell-model calculations. 
Experimentally determined strength distributions from charge-exchange reactions are needed to test 
modern shell-model calculations. We report on the measurement of the Gamow- Teller strength dis- 
tribution in ''^Co from the "'^Ni(t,''H[e) reaction with a secondary triton beam of an intensity of ~10^ 
pps at 115 MeV/nucleon and a resolution of ~250 keV. Previous measurements with the ^*Ni(n,p) 
and the ^*Ni(d,^He) reactions were inconsistent with each other. Our results support the latter. We 
also compare the results to predictions of large-scale shell model calculations using the KB3G and 
GXPFl interactions and investigate the impact of differences between the various experiments and 
theories in terms of the weak rates in the stellar environment. Finally, the systematic uncertainties 
in the normalization of the strength distribution extracted from ^*Ni(^He,t) are described and turn 
out to be non-negligible due to large interferences between the AL = 0, AS = 1 Gamow- Teller 
amplitude and the AL — 2, AS = 1 amplitude. 



PACS numbers: 26.50.-|-x,25.45.Kk,25.55.-e,27.40.+z,21.69.Cs 



I. INTRODUCTION 



Weak interactions play an important role in a variety 
of astrophysical phenomena. In particular the evolution 
of massive stars during their pre-supernova stage (type-II 
or core-collapse supernovae) is strongly affected by elec- 
tron capture and /3 decay rates. Electron-capture leads 
to deleptonization of the stellar environment. The dy- 
namics of the collapse process is modified and the iso- 
topic composition of the star is changed [1] . When the 
electron-to-baryon ratio has decreased sufRciently, and 
more neutron-rich and heavy nuclei are produced, 
decay also becomes important. 
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In order to understand the later stages of stellar evo- 
lution, the network calculations must be performed with 
accurate weak-interaction rates; pf-shell nuclei (A^40- 
65) are especially important. Fuller, Fowler and New- 
man first parameterized electron-capture rates [2-6] by 
describing the Gamow- Teller strength in an independent- 
particle model. It is well-known, however, that the 
residual interactions between valence nucleons lead to a 
quenching and fragmentation of the strength [7]. For 
pf-shell nuclei, large-scale shell-model calculations [8, 9] 
were performed to estimate Gamow- Teller strengths, 
which were then used to predict weak-interaction rates 
in the stellar environment. Use of these new rates, in- 
stead of the ones calculated by Fuller, Fowler and New- 
man, have a strong effect on the predictions for the late 
evolution stages of stars [10-12]. It is, therefore, im- 
portant that the theoretical strength distributions are 
reliable and in agreement with experimental data. In 
this paper, we test theoretical predictions of the Gamow- 
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Teller distribution in ^*Co using the ''^*Ni(t,'^Hc) reaction 
at 115 MeV/nucleon. Here, the Gamow- Teller strength 
is defined so that the strength B(GT)=3 for the decay of 
a free neutron. 

Charge-exchange reactions at intermediate energies 
(E>100 MeV/nucleon) have been used widely to probe 
the spin-isospin response of nuclei (see e.g. [13, 14]). 
Both (p,n)-type (AT^, = -1) and (n,p)-type (AT^ = +1) 
reactions have been employed. Gamow- Teller transition 
strengths (AL = 0, AS = 1) probed by nuclear reactions 
can be directly connected to weak-interaction strengths 
(for (3 decay and electron capture) since the transitions, 
mediated through the aT± operator, are between the 
same initial and final states in the two processes. 

Systematic studies of Gamow- Teller transition 
strengths were first performed at lUCF [15-19] using 
the (p,n) (AT^ = —1) reaction. It was shown [20] 
that there is a proportionality between the B(GT) 
values and the measured differential cross section at 
zero momentum transfer. This relationship between 
strength and cross section also provided the foundation 
for using other charge-exchange reactions, in particular 
the (^He,t) reaction (see e.g. [21]), to measm-c Gamow- 
Teller strengths. Presently, the most extensive ('^He,t) 
program is performed at RCNP [22] (E(3He)=140-150 
MeV/nuclcon). Using the dispersion- matching tech- 
nique, an excitation-energy resolution of ~ 30 keV has 
been routinely achieved, see e.g. Refs. [23-27]. 

In the ATx = +1 direction, the (n,p) reaction has 
been used to measure Gamow- Teller strength distribu- 
tions [28-30]. Because resolutions are limited (~1 MeV), 
alternative probes have been developed. The (d,^He) re- 
action was used at RIKEN [31], Texas A&M [32], and 
the KVI [33]. The most extensive program is performed 
at the latter institution [34, 35], where resolutions in ex- 
citation energy of ~130 keV at E(d)=85 MeV/nucleon 
are achieved. 

In a recent paper [36], we showed that the (t,^He) re- 
action at 115 MeV/nucleon can also be used to extract 
Gamow- Teller strength distributions in the AT^ = +1 
direction. A study of the ^^Mg(t,'^He) reaction was per- 
formed and combined with results from the ^^Mg(^He,t) 
reaction. The details of the reaction mechanism were 
investigated. Gamow- Teller strengths extracted from 
(t,^He) and (^IIe,t) were compared with results from 
(d,^He) and (p,n), respectively, displaying a good over- 
all correspondence. Systematic errors in the extraction 
of Gamow- Teller strength using its proportionality to 
cross section at zero momentum transfer were studied 
and quantified. It was demonstrated that such errors (on 
the level of 10-20% in the case of ^^Mg) are largely due to 
interference between the AL = 0, AS = 1 Gamow- Teller 
amplitude and the AL = 2, AS = 1 amplitude medi- 
ated mainly by the tensor-r component of the interac- 
tion. The interference can be constructive or destructive 
and it was found that the error in the strength averaged 
over many states became very small. 

On the basis of these studies, we now present the re- 



sults from a measurement of the '''*Ni(t,"^He) reaction 
and report the strength distribution in the pf-shell nu- 
cleus '"'^Co. In the past, the Gamow- Teller strength 
distribution has been obtained from ^^Ni(n,p) [30] and 
58Ni(d,2He) [37, 38] reactions. The reported Gamow- 
Tcller strength distributions were inconsistent; the inte- 
grated strengths up to excitation energies of 10 MeV were 
similar, but the details of the strength distributions were 
different. 

In the stellar environment, electron-capture rates de- 
pend strongly on details of the Gamow- Teller distribu- 
tion at low excitation energies due to phase-space effects. 
In Ref. [37] it was shown that the calculated electron- 
capture rates based on the (d,^He) and (n,p) results were 
significantly different (^^20%) for conditions in the core 
of a 25 solar-mass pre-supernova star following silicon 
depletion. A discrepancy in rates for one particular nu- 
cleus may not strongly affect the overall evolutionary 
track since typically electron-capture rates for about 5 
nuclei arc important at any stage [11]. It is important, 
however, to resolve the present ambiguity for the Gamow- 
Teller strength distribution in ^^Co, since only a limited 
number of experiments arc available to validate the shell- 
model calculations that are used to estimate the weak 
transition rates for a wide variety of nuclei. 

We note that ■''^Ni(t,^IIe) experiments have been per- 
formed at beam energies of 8 MeV/nucleon [39] and 40 
MeV/nucleon [40]. At these lower energies, systematic 
errors in the extraction of Gamow- Teller strengths can 
become large, due to the large magnitude of the tensor- 
T interaction [41] and contributions from two-step pro- 
cesses [22] , and are thus not used to validate the Gamow- 
Teller strength distributions predicted by shell models. 



II. EXPERIMENT 

The production of a secondary triton beam at NSCL is 
described in detail in Refs. [42,43]. The methods specific 
to the current data set include a primary a beam of 140 
MeV/nucleon, accelerated in the K1200 cyclotron. The 
beam impinged on a thick Be production target (9.25 
g/cm^). Tritons with a mean energy of 115 McV/nucleon 
and an energy spread of 1% were selected in the A1200 
fragment separator and guided onto a target positioned 
at the object point of the S800 magnetic spectrometer 
[44]. The spectrometer is used to detect "^He particles 
from the (t,'^He) reaction and is operated in dispersion- 
matched mode to correct for the energy spread of the 
triton beam. It was positioned at 0° degrees and the 
solid angle covered in the measurements was about 20 
msr. 

The S800 focal-plane detector [45] was composed of two 
2-dimensional cathode-readout drift detectors (CRDCs). 
They are used to measure the position and angles in the 
dispersive direction and non-dispersive direction at the 
focal plane. The positions and angles were calibrated 
by inserting masks with holes and slits at known lo- 
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cations in front of the CRDCs and illuminating them 
with a '^He beam [44]. For ray-tracing purposes, the ion- 
optical code COSY Infinity [46] was used to calculate 
the ion-optical transfer matrix of the S800 spectrome- 
ter [47] from the measured magnetic field maps. Ma- 
trix elements up to fifth order were used in the recon- 
struction of 5 = {E — Eq)/Eq {Eq is the kinetic en- 
ergy of the particle following the central-ray trajectory 
through the spectrometer and E the energy of the mea- 
sured particle). The track angles in the dispersive and 
non-dispersive directions were also obtained in the ray- 
tracing procedure and used to calculate the ^He scatter- 
ing angle. The angular resolution was 10 mrad (FWHM). 
Prom these reconstructed parameters, the excitation en- 
ergy was obtained from a missing-mass calculation with 
a resolution in the case of the ^^Co spectrum of about 
250 keV (FWHM). Angular distributions were extracted 
up to center-of-mass scattering angles of 5°. 

Two thin plastic scintillation detectors (El and E2) 
were positioned behind the CRDCs. The fast timing sig- 
nal from the El scintillator served as the trigger for the 
data acquisition, as well as the start of the time-of-fliglit 
(TOF) measurement. The stop signal for the TOF mea- 
surement was provided by the cyclotron radio frequency 
(RF). The Ught-output signal of El and E2 and the TOF 
signal were used to identify the ^He particles in the focal 
plane. 

The triton-beam intensity was monitored using an in- 

bcam scintillator (IBS) placed 30 m upstream from the 
target. The triton beam of 115 MeV/nuclcon cannot be 
bent into the focal plane of the S800 spectrometer due to 
its high rigidity. The transmission from the IBS to the 
target was, therefore, measured by comparing the rates 
measured at the focal-plane scintillator and the IBS us- 
ing a secondary ^He beam without a target and found 
to vary between 80 and 95%. The secondary beam fills 
a large fraction of the available phase-space and, as a 
result, small changes in ion-optical settings or primary 
beam properties can lead to significant changes in trans- 
mission. On average, the triton beam intensity at the 
target was about 1x10^ pps. 

Besides the measurements using a 99.93% isotopically 
enriched ^*Ni target with a thickness of 7.61 mg/cm^, 
measurements were also performed on a CH2 target with 
a thickness of 6.72 mg/cm^. The well-known Gamow- 
Teller transition to the ^^B ground state provides a con- 
venient way to check the procedures for extracting angu- 
lar distributions. Because of the relatively strong kine- 
matical correlation between momentum and angle for re- 
actions on ^^C, the excitation-energy resolutions (~ 300 
keV below 1° to ^ 450 kcV above 4°, respectively) were 
slightly worse than those obtained in the data taken on 
the 5^Ni target. 

In the experiment using the '"'^Ni target, the down- 
stream CRDC detector was partially (~15%) inefficient 
in the dispersive direction for ^He particles. Significantly 
reduced efficiency was confined to regions in the detector 
corresponding to excitation energies in ^*Co higher than 



10 MeV. Although momentum is mostly determined from 
the position measurement in the upstream CRDC (it is 
located near the true focal plane of the spectrometer), 
the scattering angle for events with missing downstream 
CRDC signal could not be well-determined. Since re- 
construction of the angular distributions is important for 
extracting Gamow- Teller strengths from the spectra, the 
detailed analysis of the data taken on the ^^Ni target was 
restricted to excitation energies below 10 MeV in ^^Co. 



III. THE ^2C(t,^He)^^B REACTION. 

In Fig. la, the excitation-energy spectrum from the 
CH2(t,^He) reaction is shown. In the missing-mass cal- 
culation, a ^^C target is assumed and the events stem- 
ming from reactions on hydrogen in the target appear 
as an asymmetric peak at negative excitation energies 
(The Q-value for i2c(t,3He)i2B,^gj.o^n(j g^^te) is -13.35 
MeV and for H(t,3He)n -0.764 MeV). Besides the ^^B 1+ 
groimd state, other known states in the -"^^B spectrum are 
indicated in Fig. la as well, along with their known J'^ 
assignments. The broad peak at around 7.5 MeV is due 
to several states, excited predominantly via dipole tran- 
sitions {J^ = 1^,2^) but containing components with 
J'^ = 3^ as well [48]. Besides the ground state, two 
other 1+ states in ^^B are known to exist (at 5.00 and 
6.6 MeV [48]) but were not identified in the current data 
set. 

In Fig. lb, a scatter plot of the laboratory scattering 
angle of the ^He versus the excitation energy is displayed. 
Because ^^C is assumed to be the target in the missing 
mass calculation, the recoil energy for reactions on hydro- 
gen in the target is underestimated, hence the remaining 
correlation between energy and angle. In the angular 
range imder consideration {Oiab < 4.5°), the events from 
hydrogen in the target do not interfere with the analysis 
of the 12 c data. 

The code DW81 [49] was used to calculate differen- 
tial cross sections for the transition to the "'^^B J'^ = 1+ 
ground state. The B(GT) for this transition equals 
0.990±0.002 [50], which is calculated from measured /3 
decay log/t values [48] . The optical-potential parameters 
for the outgoing ^He particle were obtained from elastic 
scattering of ^He on ^^C at 150 MeV/nucleon [51]. Fol- 
lowing Refs. [52, 53], the optical potential parameters 
for the triton were obtained by scaling the "^He poten- 
tial depths by 0.85 [53], while keeping the radii and dif- 
fusenesses constant. The one-body-transition-densities 
(OBTDs) were calculated using the code OXBASH [54] 
with the CKII interaction [55] in the p model space. The 
Gamow- Teller strength obtained in this calculation for 
the transition to the ground state (B(GT)=0.98) is 
close to the experimental value from /? decay. Wave func- 
tions of the target and residual nucleus were calculated 
in a Woods-Saxon potential. Single-particle binding en- 
ergies were determined in the code OXBASH [54] us- 
ing the Skyrme SK20 interaction [56]. The projectile- 
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FIG. 1: a) Excitation-energy spectrum for the (t,^He) reac- 
tion on a CH2 target, for which kinematic corrections were 
made assuming a ^'^C target. The peak at 'negative' excita- 
tion energies is due to reactions on hydrogen. The spectrum 
at Ex < —2 MeV has been down-scaled by a factor of 10. 
b) Scatter plot of the ^He laboratory scattering angle versus 
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FIG. 2: The differential cross section of the 
^^C(t,^He)^^B(g.s.) transition. The solid line is the 
result of the DWBA calculation (see text). 



uncertainty only. In the comparison between data and 
theory, the magnitude of Vo-r was determined by a fit 
to the data. Since the ratio Yar/^Tr was fixed, Vtt 
was adjusted accordingly. The fitted value of V^t was 
4.6 ±0.1 MeV [x^/Nd.o.f. = 1.3). The good correspon- 
dence between experiment and theory confirms the find- 
ings from the experiment on ^^Mg [36], namely that the 
DWBA calculations can be used to predict reliably an- 
gular distributions. For the analysis of the ^^Ni data, the 
value of V,- was set to 1.6 MeV using the known ratio of 
V^r/Vr = 2.9 [59]. 



target interaction was described by an effective '^He- 
nucleon interaction with spin-isospin (Vctt), isospin (V,-) 
and tensor- isospin (Vtt) components [53]. The former 
two are described by single Yukawa potentials, for which 
the ranges arc fixed to that of the onc-pion exchange 
potential, R = 1.414 fm. The tensor-isospin term is rep- 
resented by a potential of the form x Yukawa (with 
range R = 0.878 fm), multiplied by the tensor operator 
S12 = '"^i"'^ ''^ ~(Ji ■ (72- In principle, a spin-orbit term 
(Vl5t) should also be included, but it is usually set to 
zero since it was predicted [41] and confirmed by experi- 
ment [52, 57] that it hardly contributes at low momentum 
transfers. Since spin-transfer is required, only V^t and 
Vtt arc important for Gamow- Teller transitions. The 
ratio Vctt/Vtt was fixed to 1.0 following Ref. [58]. 

In Fig. 2, the differential cross section for the tran- 
sition to the '^^B ground state is compared with the 
DWBA calculations. The error bars reflect statistical 



IV. THE ^®Ni(t,^He)^*Co REACTION. 

In Fig. 3a, the ^*Ni(t,^He) spectrum is shown inte- 
grated over the full angular range covered in the present 
experiment. A binning of 250 keV is applied, correspond- 
ing to the energy resolution. The Q-value for the tran- 
sition to the J'^ = 2+ ground state of ^^Co is -0.363 
MeV, less negative than that of the H(t,'^He)n reaction. 
By examining the correlation between scattering angle 
and ^He energy, no significant contributions from events 
due to the H(t,^He)n reaction to the ^^Co spectrum were 
foimd. 

Transitions of various multipolarities contribute to the 
spectrum shown in Fig. 3a. Whereas transitions with 

AL = are associated with differential cross sections 
that are strongly forward-peaked, differential cross sec- 
tions of dipole transitions {AL = 1) peak at about 3° 
and transitions due to higher angular-momentum trans- 
fer are relatively flat at very forward angles and have a 
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FIG. 3: a) The ^*Co excitation -energy measured via the 
^*Ni(t,^Hc) reaction at 115 McV/nuclcon integrated over the 
full measured solid angle, b) The same reaction, but gated 
on events with diab =0°-l°. Tlie two indicated peaks can be 
identified as Gamow- Teller transitions because of their max- 
ima at forward scattering angles, c) The same reaction, but 
gated on events with Oiab =2°-3°. The indicated broad struc- 
ture broad is likely due to dipole transitions, identified by its 
maximum in this angular range. 



maximum at ~ 4°. In Fig. 4 angular distributions calcu- 
lated in DWBA for several relevant transitions to states 
of different multipolarity aic compared. The details of 
these calculations are discussed below. 

To illustrate the method of using the angular distribu- 
tions to identify different types of transitions in the data, 
the excitation energy spectra for 0° < 0iab(^He) < 1° 
and 2° < 6'iab(^He) < 3° are displayed in Figs. 3b and 
3c, respectively. Just below = 2 MeV and at about 
Ex = 4, MeV two strongly forward-peaked transitions can 
be identified (as indicated in Fig. 3b), reveahng the pres- 
ence of transitions with AL — 0. For 7 < E^ < McV, 
the differential cross section rises at backward angles, re- 
vealing the presence of dipole contributions presumably 
due to the various components of the isovector spin giant 
dipole resonance (AL = 1, AS' = 1, J'' = 0-,l-,2-) 
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FIG. 4: Differential cross sections calculated in DWBA for 
^*Ni(t,^He) transitions to states in ^®Co with different J". In 

all cases, an excitation energy of 2 MeV in "'^Co is assumed 
and the normalizations are adjusted so that the cross sections 
are of the same order of magnitude for easy comparison. 



and its non-spin-flip partner, the isovector giant dipole 
resonance (AL = 1, A5 = 0, J'^ = The main bump 
is indicated in Fig. 3c. Note that even though for about 
15% of the data above E^ > 10 MeV the angle was not 
well determined in the analysis (see Section II) such gross 
features remain visible. 

Since transitions with different angular momentum 
transfer are not necessarily separated in energy, a multi- 
pole decomposition analysis (MDA) must be performed 
to extract the strength distributions as discussed below. 
In the Tz = +1 direction, the possible AL =- transi- 
tions arc to Gamow- Teller states and the Ihjj isovector 
(spin-flip) giant monopole resonances [82] . The latter res- 
onances only contribute to the '"^^Co spectrum at E^ ^ 10 
MeV [40, 60, 61] and peak at about E^ = 25 MeV [83]. 
Therefore, extraction of the AL = component in the 
region E^ < 10 McV based on angular distributions is a 
very selective tool for identifying the contributions due 
to Gamow- Teller transitions. 

A MDA was performed for each 250-keV wide bin, by 
fitting the experimental differential cross section to a lin- 
ear superposition of theoretically predicted angular dis- 
tributions with the code DW81 [49] for various types of 
transitions. The optical potential parameters for the ^^Ni 
target were taken from fitting ^He elastic scattering data 
at 150 MeV/nucleon [51]. Wave functions projected on a 
complete Ip-lh basis were calculated in a normal- modes 
procedure [62, 63] using the code NORMOD [64]. In this 
formalism, the full strength (i.e. 100% of the non-energy- 
weighted sum rule) associated with the single-particle 
multipole operators is exhausted. The calculated cross 
sections, therefore, do not relate to any particular state 
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in the spectrum and the main purpose of the calcula- 
tions was only to determine the shape of the angular 
distribution used in the MDA. The DWBA calculations 
were performed for each excitation-energy bin separately 
to account for the increase in momentum transfer with 
increasing excitation energy at a fixed scattering angle. 

The statistical accuracy of the data was limited, and 
only five l°-wide angular bins between 0° and 5° in the 
center-of-mass were used. With a limited number of 
data points, a fit with multiple components is difficult. 
Therefore, we used only two multipole components per 
excitation-energy bin: a Gamow- Teller component and 
one component with larger angular-momentum transfer. 
The choice for the second component was made based 
on the excitation-energy of the bin. In the region be- 
low E^^b MeV, states with J"" = 1+, 2+, 3+, 4+, 5+ are 
known to reside and contributions from dipolc transitions 
are small [65] . At beam energies above 100 MeV/nucleon 
and at forward angles, transitions with large orbital an- 
gular momentum transfers arc strongly suppressed and 
states with J'^ = 4+ , 5+ are hardly excited and can safely 
be ignored [36]. The transitions to the J'^ = 2+, 3+ states 
have similar distributions (sec Fig. 4). As a result, differ- 
ences between fits with a Gamow- Teller component and 
either a J'^ = 2+ or a J'^ = 3+ component resulted in 
difi^erences for the extracted Gamow- Teller cross sections 
that were smaller than the statistical error. Therefore, 
for extracting Gamow- Teller strength in the excitation- 
energy region below 5 MeV, the MDA was performed 
using a Gamow- Teller and J'^ = 2+ component only. As 
an example of the procedure, the result of the MDA in 
the energy bin from 1.75 to 2 MeV is shown in Fig. 5. 
A systematic error of 5% was assigned to the extracted 
Gamow- Teller strength in each bin based on the differ- 
ences between fits with a J"^ = 2+ or a J'^ = 3+ compo- 
nent. 

For the region 5 < Ex < 10 MeV, the situation is 
more complicated, since dipole transitions will play an 
increasingly larger role at higher excitation energies. It 
was found that the fits to the angular distributions had 
consistently lower values per degree of freedom when 
using a dipole component [84] in the MDA in addition 
to a Gamow- Teller component, compared to using tran- 
sitions leading to states with = 2+ or = 3+. Irre- 
spective of the choice for the component used in addition 
to the contribution from the Gamow- Teller component, 
the extracted Gamow- Teller cross section did not change 
significantly (i.e. by more than statistical error bars). A 
systematic error of 10% of the extracted Gamow- Teller 
in each bin was assigned based on the differences between 
fits with a J'^ = 2+, 3+ or a dipole component. 



^*Ni(t,^He) 

E(t)=115 MeV/nucleon 
1 .75 MeV<E^{^^Co)<2.0 MeV 

data 
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FIG. 5: The differential cross-sections for the energy bin 
1.75 < E^i^^Co) < 2.0 MeV and the result of the MDA (solid 

line) using a linear combination of 1+ (dashed line) and 2'^ 
(dotted line) components. The error bars in the data are of 
statistical nature only. 



energy spectrum to strength. This is done by using the 
proportionality between the cross section at zero momen- 
tum transfer {^{q = 0)) and Gamow- Teller strength 
(B(GT)) as derived in Eikonal approximation [16, 20]: 



da 

In 



{q = Q)=KN''\J,r?B{GT). 



(1) 



Here, K 



where £^i(/) is the reduced energy 

in the incoming (outgoing) channel. | Jo-t|^ is the volume 

integral of the (tt component of the projectile-target in- 
teraction. The distortion factor is defined by [20] 



N 



dnpWBA 



(2) 



and is determined using the DWBA code. For the plane- 
wave (PWBA) calculation, the depth of the optical po- 
tentials and the charges of the target and residual nuclei 
are set to zero. The factor KN^\J„t\^ is referred to as 
the unit cross section a. 

Apphcation of Eq. 1 to calculate B(GT), requires the 
knowledge of the experimental cross section at zero mo- 
mentum transfer (i.e. requiring that both the Q-value 
(Q) of the transition and the scattering angle are zero), 
which is obtained from: 



EXTRACTION OF GAMOW-TELLER 
STRENGTHS. 



da 

dn 



iq = 0) 



ln(g = o) 
n(Q'0°) 



da 



DWBA 



da 

dn 



iQ,o°) 



(3) 



The next step in the analysis is to convert the ex- 
tracted Gamow- Teller contributions in the excitation- 



In this equation, 'DWBA' refers to values calculated in 
the DWBA code. The experimental cross section a,t 9 = 
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0° and its error are taken from the fitted Gamow- Teller 
angular distribution in the MDA. 

In practice, unit cross sections are determined using a 
transition for which B(GT) is known from (3 decay; we 
apply the same method here, albeit indirectly. As with 
the procedure used in the analysis of the ^^Ni(d,^He)^^Co 
experiment [37, 38] the unit cross section was determined 
from the transition to the strongest 1+ state at = 1.87 
MeV present in the spectrum. The analog of this transi- 
tion has been studied with the ^*Ni('^He,t)^^Cu reaction 
at 140 MeV/nucleon [66, 67]. Under the assumption of 
isospin symmetry, the difference in B(GT) between the 
analogs is only the square of an isospin Clebsch-Gordan 
coefficient. In the case of the ^^Ni target, the strength is 
a factor of 6 larger for the AT^ = +1 transition than for 
the ATz = — 1 transition. The Gamow- Teller strength for 
the ground state transition from ^^Ni to ^*Cu was deter- 
mined from the measured log ft value for P decay [65]. 
This yields a Gamow- Teller strength of the l"*" state at 
1.87 MeV in ^SCo of 0.72 ± 0.05. Another 1+ state is 
known to be located at = 1.73 MeV in ^^Co [65]. Its 
B(GT) was determined to be 0.17±0.04 [67]. Since in the 
present (t,^He) experiment this state was not separated 
from the one at E^ = 1.87 MeV, the unit cross section 
was determined using the sum of these states. This cal- 
ibration procedure has a systematic error that will be 
discussed in more detail in Section VII. Once the unit 
cross section was determined, it was used to convert the 
cross sections at 0° obtained from the MDA analysis in 
each energy bin to strengths, using Eqs. 1-3. 



VI. RESULTS AND COMPARISON WITH 
THEORY AND PREVIOUS DATA. 

In Fig. 6 the extracted Gamow- Teller distribution for 
excitation energies between MeV and 10 MeV in ^*Co 
from the ^*Ni(t,'^IIe) experiment is shown and compared 
with theory. In both theoretical calculations, a quenching 
factor of (0.74)^ [68] is appHed. The different interactions 
result in significantly different strength distributions. 

To facilitate the comparison between the data and the- 
ory in Fig. 6b, the calculations shown in Fig. 6a were 
folded with the experimental resolution and binned in 
the same manner as the data. The calculation with 
the KB3G interaction does relatively well in describing 
qualitatively the experimental strength distribution up 
to about Ex = 4: MeV, but is too large in magnitude. 
Above Ex — 4 MeV, too little strength is predicted. The 
calculation with the GXPFl interaction does not repro- 
duce the strongest state at 1.87 MeV, but predicts more 
strength at higher excitation energies and is in that re- 
spect more consistent with the data. 

Next, we compare the results from the ^*Ni(t,'^IIe) 
experiment with previous experimental results from 
58Ni(n,p) [30] and 58Ni(d,2He) [37, 38]. Better energy 
resolution is achieved in the latter experiment (130 keV). 
The energy resolution in the (n,p) experiment was 1.3 
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SM GXPF1 
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FIG. 6: a) Large-scale shell-model calculations for the 
Gamow- Teller strength distribution in ^®Co using the KB3G 
[69] interaction [37, 38, 70] (black) and GXPFl [71, 72] in- 
teraction [38] (green), b) Gamow- Teller strength distribution 
extracted from the ^*Ni(t,"^He) data compared with the the- 
oretical distributions that were folded with the experimental 
resolution (250 keV) and binned in the same manner as the 
data. 



MeV. Therefore, when comparing the various results, the 
resolution must be taken into account, as in the compar- 
ison between the '^^Ni(n,p) and ^^Ni(d,^IIe) results made 
in Refs. [37, 38]. However, in those Refs. the bins for the 
two data sets were shifted by 0.5 MeV relative to each 
other, which made a direct comparison somewhat difh- 
cult. Therefore, we performed the comparison between 
the (d,^He) and (t,'^He) results and the comparison be- 
tween the (n,p) and (t,'^He) results separately. 

In Fig. 7a, the (d,^IIe) and (t,''He) results are com- 
pared. Results from ^^Ni(d,^IIe) were directly taken from 
Table 2 in Ref. [37]. For the 5*Ni(t,3He) data, two al- 
ternative methods to extract the Gamow- Teller strength 
distribution in 1-MeV bins were used: i) The MDA analy- 
sis was performed in each 1-MeV bin and ii) the extracted 
B(GT)s from four consecutive 250-keV bins were added. 
The differences between methods i) and ii) were smaller 
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than the statistical errors, and in Fig. 7a the results from 
method ii) are shown. Since the energy resolution in both 
data sets is much smaller than the bin size, the difference 
in energy resolution hardly affects the representation and 
was not corrected for. The (d,^He) and (t,^He) data sets 
are consistent. The x^/A^d.o.f. — l-2f with A^d.o.f = 9, i.e. 
well within a 95% confidence interval (since the absolute 
scale of the distribution was determined in the same man- 
ner for the two data sets, A^d.o.f = 9 instead of 10). Also 
included in the figure are the theoretical distributions, 
folded with the experimental resolution achieved in the 
(t,'^He) experiment and binned in the same manner as 
the data. 



In Fig. 7b, the (n,p) and (t,'^He) results are compared. 
Again, a bin size of 1 MeV is used, but shifted by 0.5 
MeV relative to Fig. 7a. Since the (n,p) data had a 
much lower energy resolution, the resolution of (t,'^He) 
data was artificially reduced to that of the (n,p) data. 
The two different theoretical calculations were also folded 
with the resolution from the (n,p) experiment and binned 
according to the data. As can be seen from the (t,''He) 
strength distribution in Fig. 7b, the folding and choice 
of binning almost hides the presence of a strong peak at 
Ex — 1.87 MeV, and it is, therefore, understandable that 
it was not seen in the (n,p) data. However, the (t,'^He) 
and (n,p) data sets are not consistent {x^/N^.o.i. = 6.0 
with iVd.o.f = 10, i.e. well outside the 99.5% confidence 
interval). The unit cross section in the (n,p) experiment 
[30] was obtained by studying its systematic behavior in 
nuclei of similar mass. Although this could lead to a sys- 
tematic discrepancy with the (t,'^He) results, we note that 
rescaling either of the data sets does not improve the con- 
sistency, since for < b MeV (E^ > 5 MeV) the (n,p) 
strengths are larger (smaller) than the (t,'^He) strengths. 
It is important to note that, if the data sets with rela- 
tively good resolutions (from (d,^IIe) and (t,'^IIe)) were 
not available, one might have concluded that the theo- 
retical calculation with the GXPFl interaction is favored 
for predicting the Gamow- Teller distribution since, ex- 
cept for a minor shift in excitation energy, the match 
with the (n,p) results is good. But after comparison with 
the (d,'^IIe) and (t,'^IIe) results it becomes clear that the 
strong state at E^ = 1.87 MeV is missed when using the 
GXPFl interaction. Reiterating a conclusion from Ref. 
[37], it is necessary to obtain data with good resolution 
to test shell-model calculations in detail. 
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FIG. 7: a) Comparison of the results of the ^**Ni(d,^He) and 
^*Ni(t,''He) experiments and the theoretical predictions. A 
binning of 1 MeV was applied and the theory was folded 
with the experimental resolution of the (t,''He) experiment 
(250 keV) before binning, b) Comparison of the results of 
the ^*Ni(n,p) and ^*Ni(t,'^He) experiments and the theoreti- 
cal predictions. A binning of 1 MeV was applied. Note the 
0.5 MeV shift relative to a). The (t,^He) data set and the- 
ory were folded with the experimental resolution of the (n,p) 
experiment (1.3 MeV) before binning. 



VII. THE NORMALIZATION OF B(GT) 



Although the two theoretical calculations and the 
three data sets have different strength distributions, the 
summed strengths up to E^ — 10 MeV are very sim- 
ilar: Ekb3G [B(GT)] = 4.0, Egxpfi [B(GT)] = 4.1 
and E(t.3Ho) [B(GT)] = 4.1 ± 0.3, E(d,^He) [B(GT)] - 
3.4 ± 0.2 and J^in,?) [B(GT)] = 4.1 ± 0.1 (the errors in 
the summed strengths for the three data sets are statis- 
tical only). 



As discussed in Section V, the Gamow- Teller strength 
distribution in ^^Co was normalized using known 
strengths for the analog of the 1+ states at 1.73 MeV 
and 1.87 MeV studied with the ^*Ni(3He,t) reaction at 
140 MeV/nucleon [66, 67]. In turn, the B(GT) of these 
analog states, located at 10.60 MeV and 10.82 MeV in 
^*Cu, were calibrated with the B(GT) known from /3 de- 
cay of the ground state transition ''^^Ni(g.s.)-^^^Cu(g.s.). 
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Transitions from '''^Ni to '"^^Cu have also been studied 
with the (p,n) reaction [18] at 120 MeV and 160 MeV. 
After determining the unit cross section for the ground- 
state transition, in both the (^He,t) and (p,n) data sets 
the B(GT) for the second 1+ state at 1.05 MeV was ex- 
tracted as well. However, the ratios of B(GT)s for the 
transition to the ground state and to the second 1+ state 
in ^*Cu are very different for the three data sets. This 
is shown in Table I. In the (p,n) experiments, the '''^Cu 
ground state was not separated from the isobaric analog 
state {J'^ = 0+) at ^^(^^Cu)=0.20 MeV, which could 
have led to a systematic error. Perhaps, the difference 
between the ratios between the ratios in Table I for the 
two (p,n) experiments is an indication of that. However, 
this cannot fully explain the large difference between the 
two (p,n) and the (^He,t) result. 

In Ref. [73], the Gamow- Teller strength distributions 
in 58 Cu calculated using the KB3G and GXPFl inter- 
actions were discussed. When using the GXPFl inter- 
action, the overall strength distribution was well repro- 
duced, but the B(GT) for the transition to the ground 1+ 
state at Ex=1.05 MeV was too low compared to exper- 
iment. With the KB3G interaction, the B(GT) for that 
state was better reproduced, but a group of 1+ states 
experimentally found around £^=3.5 MeV was absent. 

In a simple independent-particle shell model, the neu- 
trons and protons in ^^Ni fill all orbitals up to /7/2- The 
two remaining neutrons populate the P3/2 and ^5/2 or- 
bitals. Therefore, relatively strong contributions from 
^'P3/2-^./5/2 and J^/5/2-7i'P3/2 particle-hole components 
are to be expected for excitations of the lowest lying 1+ 
states in "''^Cu, even though these components arc purely 
AL = 2, AS = 1 nature (i.e. not Gamow- Teller). As 
discussed in Ref. [73], in more realistic models using 
the GXPFl and KB3G interactions, this is still largely 
true: transitions to the 1+ ground state and excited state 
at Ex~l.05 MeV have strong contributions from the p- 
orbit. The components with AL = 2, AS = 1 tend to 
break the proportionality between cross section at zero 
momentum transfer and B(GT) in Eq. 1 [20, 36]. Since 
the effect of the interference between AL = 2, AS = 1 
and AL = 0, AS = 1 amplitudes on the cross section de- 
pends on the reaction and beam energy, it could possibly 
explain the discrepancies between the (■^He,t) and (p,n) 
results shown in Table I. 

In the i.e. (n,p) direction, contributions from 
AL = 2, AS = 1 components to the Gamow- Teller tran- 
sitions are relatively small in the case of the '"'^^Ni target, 
since (in an independent-particle picture) the protons do 
not fill the P3/2 and /5/2 orbitals. 

To estimate the proportionality-breaking effects we 
performed cross-section calculations for the 5^Ni(p,n) 
and 5*Ni(''Hc,t) reactions for the transitions to the 1+ 
ground state and excited state at E^ — l-OB MeV in ^^Cu. 
The beam energy was set to 140 MeV/nucleon, which 
is between the beam energies of the two available (p,n) 
measurements and equal to the energy used in the (^He,t) 
experiment . 



For the 5^Ni(p,n) reaction, the code DW81 [49] was 
used. The Lovc-Francy effective interaction at 140 MeV 
[41, 74] was employed and exchange effects were treated 
exactly. OBTDs for the transitions to the relevant 1+ 
states in ^^Cu were taken from shell-model calculations 
in which at most 5 nucleons are allowed to be excited 
from the /5/2 orbit to higher-lying orbits [73, 75]. To test 
the dependence of the results on the interaction, two sets 
of OBTDs were used, calculated with the KB3G and the 
GXPFl interactions. Radial wave functions of the target 
and residual nuclei were calculated using a Woods-Saxon 
potential. Single-particle binding energies of the parti- 
cles were determined in OXBASH [54] using the Skyrme 
SK20 interaction [56]. Optical-model parameters were 
calculated following the procedure in Ref. [76] which in- 
cludes Coulomb corrected isovector terms to account for 
the differences between the incoming (proton plus ^^Ni) 
and outgoing (neutron plus ^^Cu) channel. 

In order to perform a precise comparison with the (p,n) 
results, the calculations for the 5^Ni(^He,t) reaction need 
to be performed with the same effective interaction (i.e. 
the Love-Franey interaction [41, 74] instead of the effec- 
tive "^Hc-nucleon interaction used above). Therefore, the 
code FOLD [77] was employed (for more details see Ref. 
[36]), in which the Love-Franey interaction is double- 
folded over the transition densities. Exchange is treated 
in the short-range approximation described in Ref. [41]. 
OBTDs for the target-residual nucleus system arc the 
same as those used in the (p,n) calculation [75] . For ^He 
and ^H, densities were obtained from Variational Monte- 
Carlo results [78]. Optical- model parameters for the in- 
coming ^ He channel were taken from Ref. [51]. Following 
the analysis in Ref. [53], the potential- well depths were 
scaled with a factor 0.85 for the outgoing triton channel. 

All cross-section calculations were performed at zero- 
momentum transfer by setting the Q-value of the reac- 
tion to zero and using the cross section at 0° scattering 
angle. For each transition and for each set of OBTDs, 
three calculations were performed: i) a plane-wave (PW) 
calculation, ii) a distorted-wave (DW) calculation and 
iii) a distorted-wave calculation in which the tensor-r 
amplitudes of the Love-Franey interaction were set to 
zero (DWtt=o)- For each calculation, a unit cross sec- 
tion was calculated: a = ^{q = 0)/B(GT) and subse- 
quently the ratio of the unit cross sections for the transi- 
tions to the ground state and the state at Ea;=1.05 MeV 
ijf^)- If this ratio equals unity, the proportionality be- 
tween B(GT) and cross section at zero degrees is perfectly 
maintained and a deviation from unity signifies propor- 
tionality breaking between the two transitions. In ad- 
dition, the ratio of the distorted-wave calculations with 
and without the tensor-r component of the effective in- 
teraction was calculated, so that it can be determined 
to what extent the breaking of the proportionality stems 
from this source. 

The results of the calculations are summarized in Table 
II. In the absence of distortions (plane- wave calculation; 
PW) the proportionality is perfect because all compo- 
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TABLE I: Comparison between Gamow- Teller strengths extracted for the ground state and second 1"*" state in ^*Cu via ^*Ni(p,n) 
at 120 MeV and 160 MeV [18] and ^®Ni(^He,t) [66, 67). 





^*Ni(p,n) 120 MeV " 


^*Ni(p,n) 160 MeV" 


5*Ni(3He,t) 140 MeV/nucleon 


B(GT) ""Cn (1+ g.s.) 


0.165" 


0.165" 


0.155±0.001" 


B(GT) ■'''Cu (1+ 1.05 MeV) 


0.5 


0.4 


0.265±0.013 


B(GT)'=«Cu(l+g.s.) 
B(GT)SSCu(l+1.06 MeV) 


0.33 


0.41 


0.58 



"No explicit error bars for B(GT) values are given in Ref. [18]. 
The authors note that the error in the absolute cross sections for 
each state are ~ 15% 

'Fixed to the known B(GT) from /3 decay [65]. The small dif- 
ference between the value for the (p,n) reactions and the ("^Hejt) 
reaction stems from the constants used when converting the log/t 
value to B(GT). 



TABLE II: Results of a theoretical study of the discrepancy between the ratio of the extracted B(GT)s for the transitions to 
the ground state and first-excited state at 1.05 MeV in ^**Co for the ^*Ni(p,n) and ■■''**Ni(-^He,t) reactions. The OBTDs that 
are used in the calculations axe from Ref. [75]. 'PW and 'DW refer to the cross section calculated in the plane- wave and 
distorted-wave approximation, respectively. 'DWtt=o refers to cross sections from the distorted-wave calculation in which the 
tensor-T amplitudes in the Love-Franey interaction have been set to 0. For further details, see text. 









CtXITT" 










KB:iG" 










5«Ni(p,n)'^*Cu 




68 


Ni(3He,t)-'«Cu 




^«Ni(p,n)5«Cu 




68 


Ni(='He,t)58 


Cu 




g.s. 


1.05 MeV 


CT1 .n.=i 


g.s. 


1.05 MeV 


CTi .n.=i 


g.s. 


1.05 MeV 




g.s. 


1.05 MeV 


CT1 .n.=i 


B(GT)th 


0.275 


0.318 




0.275 


0.318 




0.496 


0.768 




0.496 


0.768 




ihil = 0) 
(mb/sr) 

PW 


4.52 


5.28 


1.00 


36.9 


42.9 


1.00 


8.19 


12.8 


1.00 


66.5 


103. 


1.00 


DW 


1.61 


1.58 


1.19 


4.01 


3.24 


1.44 


2.77 


3.73 


1.15 


6.73 


7.80 


1.34 


DWTr^O 


1.21 


1.49 


0.94 


3.33 


3.62 


1.07 


2.25 


3.63 


0.96 


5.78 


8.60 


1.04 


DW/DWTr = 


1.33 


1.06 




1.20 


0.89 




1.23 


1.03 




1.16 


0.91 





"In the calculations shown in this table, the strengths and cross 
sections have not been adjusted for the Gamow- Teller quenching 
factor of 0.74^ [68] 



nents with AL = 2, AS = 1 are zero at q=0. In the full 
distorted- wave (DW) calculation with the GXPFl inter- 
action, the ratios of unit cross sections (Z^-^) are equal 

^0-1.05' ^ 

to 1.19 and 1.44 for the (p,n) and (^He,t) reactions, re- 
spectively. Those numbers are similar when using KB3G: 
1.15 and 1.34. In brief, significant proportionality break- 
ing effects are seen in both reactions, and the effect on 
the ratio of the cross sections for the excitation of the 
first two 1+ states is about twice as large for the (^He,t) 
reaction as for (p,n). Therefore, according to the calcu- 
lations with GXPFl: 



The ratio 



J(p,n) 



3(He,t) 



GXPFl 



and with KB3G: 



J(p,n) 



3(He,t) 



KB3G 



1.19 

L44 



1.15 
L34 



0.826 



: 0.858 



(4) 



(5) 



calculated from the ex- 



perimental results in Table I equals 0.33/0.58=0.57 if the 
(p,n) data taken at 120 MeV are used and 0.41/0.58=0.7 
if the (p,n) data taken at 160 MeV/nucleon are used. In 
short, the discrepancy between the (p,n) and (^He,t) data 
are qualitatively explained. It is hard to reach stronger 
conclusions about the quantitative agreement, owing to 
the 15% error bars in the (p,n) data. New ^^Ni(p,n) 
experiments at Ep = 200 MeV and 300 MeV have re- 
cently been performed [79], but not yet fully analyzed. 
The results could shed further light on the analysis. 

From Table II it can be seen that the changes in the 
unit cross section solely due to the inclusion of tensor- 
r components in the interaction are predicted to be 
slightly stronger for (p,n) than for (■'He,t) since the ra- 
tio DW/DWtt=o is higher for the transitions with the 
former probe. However, for the (p,n) reaction, these con- 
tributions counteract the proportionality breaking due 
to other causes (such as exchange and AL = 2, AS = 1 
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contributions mediated through the err component of the 
interaction) on the level of 4-6%, whereas for the ('^He,t) 
reaction they reinforce such effects. In addition, the sign 
of the interference is the same for both transitions in the 
case of the (p,n) reaction and opposite in the case of the 
('^He,t) reaction. Therefore, the ratio of unit cross sec- 
tions is more strongly affected in the case of the (•^He,t) 
reaction. 

Besides the calculations in Table II, we also checked 
the effect of removing all contributions from the vpzli' 
7r/5/2 and J^/5/2-7rp3/2 components in the cross section 
calculations. The value of iZ^-i- — 1) became 0.91 for 
the (p,n) reaction and 0.95 for the (■^IIe,t) reaction, with 
very minor differences between the results with KB3G or 
GXPFl. This confirms that these particle- hole compo- 
nents are indeed the leading cause for the discrepancies 
between the two reactions and the breaking of the pro- 
portionality. 

A systematic error in the absolute scale of the Gamow- 
Teller strengths extracted from ^*Ni('^IIe,t) directly 
translates into a systematic error in the absolute scale of 
the strengths extracted from ^*Ni(t,'^He). Since, accord- 
ing to the calculations, the cross section of the transition 
to the ground state in ^^Cu is increased due to the inter- 
ference between AL = 2, AS* = 1 and AL = 0, AS" = 1 
amplitudes by about 20-25% (DW/DWTr=o=l-2 and an 
additional 5% is included to account for other proportion- 
ality breaking effects), the B(GT) for the other states in 
the ^^Cu spectrum (including the T> used for calibrat- 
ing the B(GT) in ^*Co) are underestimated by the same 
percentage. The same effects occur for the results from 
the (dj^He) experiment for which a similar calibration 
procedure was used. 




FIG. 8; a) Electron-capture rates on ^*Ni in the stellar en- 
vironment at pYe — 10^ gcm""^ using various theoretical 
(green lines) and experimental (black lines) B(GT) distribu- 
tions in ^*Co, as labeled in c). b) Same a a) but calculated at 
pYe — 10® gcm~^. c) The ratio of electron capture rates in a) 
to the FFN rates, d) The ratio of electron capture rates in b) 
to the FFN rates. Note the difference in scales on the Y-axis 
for the plots. For details, see text. 



VIII. ELECTRON-CAPTURE RATES 

To understand how the differences between the various 
measured and calculated Gamow- Teller strength distri- 
butions affect the electron capture rates, we calculated 
these rates for various stages during stellar evolution. 
The method to calculate electron-capture rates is de- 
scribed in detail in Refs. [2-5] and implemented in a new 
code [80] . Calculations were performed in a grid spanning 
pYe values from 10^ gcm~^ to lO^'' gcm~^ and T values 
from 0.01 X 10^ K to 100 x 10^ K. Here, we present the re- 
sults for two different pY^ regimes of relevance in the late 
stages of evolution of massive stars (11-40 solar masses) 
[11]: pYe = 10^ gcm~'^, which corresponds to conditions 
during Silicon burning and Silicon depletion (Yg ~ 0.47; 
this regime was also investigated in Refs. [37, 38]) and 
pYf. — 10^ gcm""^, which corresponds to conditions dur- 
ing the pre-supernova stage (Yg ~ 0.44; ~ 0.5 s before 
core bounce). In the former case, the temperature is 
T ~ 4 X lO^K and in the latter case, T ~ 8 x lO^K. In the 
present calculations, we only consider transitions from 
the parent ground state to daughter states described by 
Gamow- Teller strength distributions and ignore transi- 



tions from thermally populated parent states. Especially 
at the higher temperatures, this leads to an underestima- 
tion of the electron-capture rates [9]. However, here we 
are mostly interested in relative deviations in the rates 
due to the differences in strength distributions. 

In addition to the original results by Fuller, Fowler 
and Newman (FFN) [2-6], electron-capture rates were 
calculated for five Gamow- Teller strength distributions 
in ^^Co: i) The theoretical prediction employing the 
GXPFl interaction [38], ii) the theoretical prediction em- 
ploying the KB3G interaction [37, 38, 70], iii) the distri- 
bution extracted from the ^*Ni(t,'^He) experiment, iv) 
the distribution extracted from the ^^Ni(d,^He) exper- 
iment [37, 38] and v) the distribution extracted from 
the ^*Ni(n,p) experiment [30]. For cases i) and ii) the 
I''" states in ^^Go were positioned at exactly their cal- 
culated values. For case iii), the strength extracted for 
each 250-keV wide bin was placed at the center of that 
bin. For case iv), the strength was distributed according 
to the values extracted state-by-state in Table 1 of Ref. 
[37] below Ej,(5®Co)=4. MeV. Above that energy, val- 
ues extracted per 1-MeV bin (Table 2 of Ref. [37]) were 
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equally distributed over the respective bins, and placed 
at the center of four 250-keV bins. Finally, for case v) 
the strengths extracted per 1-MeV bin were equally dis- 
tributed over the respective bins and placed at the center 
of four 250-keV bins. 

The results of the calculations are plotted in Fig. 8. 
Fig. 8a shows the rates calculated at the lower densi- 
ties, at temperatures T= 2 - 10 x lO^K and Fig. 8c 
displays the same calculations, but relative to the FFN 
results (note that the quenching factor of 0.74^ is not in- 
cluded in the FFN values). At this low density, and cor- 
respondingly low electron chemical potential {(j)e « 0.7 
MeV at T= 1.0 x 10^), the rates rise rapidly with tem- 
perature due to the fact that electrons can have ener- 
gies ~ ksT larger than {ksT = 90 keV (900 keV) at 
T= 1 X lO^K (10 X lO^K), with ks the Boltzmann con- 
stant) and thus increasingly populate the lowest-lying l"*" 
states. This thermal broadening of the electron Fermi 
surface is characterized by the degeneracy parameter 
(l>e/{kBT) (here defined positive) which is reduced ("lift- 
ing" of degeneracy) from « 7.7 at T= 1 x lO^K to 
« —0.35 at T= 10 X lO^K. The temperature dependences 
of the rates at the low density are, therefore, very sen- 
sitive to the precise location of these low-lying states, 
and when calculated from experimental strength distri- 
butions, also to the binning and resolution of the data. 
The rate calculation using the strength distribution ob- 
tained with the GXPFl interaction result in much lower 
rates than the one employing KB3G at this density, ow- 
ing to the near absence of strength in the capture win- 
dow (« {(j)e + rrieC^ — w) where w is the ground-state- 
to-ground-state capture threshold and nieC^ the electron 
rest mass). Except for the lowest temperatures, the rates 
calculated with the strength distribution from the (n,p) 
experiments are relatively low compared to those from 
(d,^He) and (t,^He) experiments. 

At the higher density, 0e ~ 4.7 MeV. As a result, a 
larger fraction of the strength distribution can be ac- 
cessed compared to the case at lower densities. The de- 
generacy (f>e/(kBT) reduces from «52 at T= 1 x lO^K 
to «4.66 at T= 10 x lO^K, but the Fermi-Dirac distri- 
bution still resembles a sharp energy filter that accesses 
daughter states almost exclusively in the capture win- 
dow. The dependence of rate on temperature is thus rel- 
atively weak, as shown in Fig. 8b. Rate estimates using 
strength distributions that have relatively little strength 
located within the capture window (i.e. the one extracted 
from the (n,p) experiment and the theoretical calculation 
using the GXPFl interaction) will result in lower rates 
than those exhibiting more Gamow- Teller strength at low 
excitation energies (i.e. the (t,^He) and (d,^He) exper- 
imental results and the theoretical calculation with the 
KB3G interaction), especially since the phase space avail- 
able for capture depends strongly (power of ~ 5) on the 
total electron energy. The FFN rates have a significantly 
steeper dependence on temperature. 



IX. CONCLUSIONS 

The Gamow- Teller strength distribution in ^^Co has 

been extracted from the ■''*^Ni(t,'^He) reaction at 115 
MeV/nucleon. Although the statistical errors are rela- 
tively large, the results are important to cross check exist- 
ing results from ■''®Ni(d,^He) and ^®Ni(n,p) experiments, 
which were inconsistent. Our data are consistent with 
those from the ^^Ni(d,^He) reaction and deviate from 
the ^SNi(n,p) data. 

Comparisons between the experimentally extracted 
strength distributions and shell-model calculations using 
the GXPFl and KB3G interactions were made. The pre- 
dictions with the KB3G interaction describe the strength 
distributions at excitation energies below 4 MeV in ^*Co 
well, but the calculations with the GXPFl interaction 
better reproduce the strength distribution at higher ex- 
citation energies. 

Systematic uncertainties in the calibration of the abso- 
lute Gamow- Teller strength scale were investigated. Such 
uncertainties are due to potentially large interference ef- 
fects between AL = 2, AS' = 1 and AL = 0, AS = 
1 components for the transition ^^Ni(''Hc.t)'^^Cu(g.s.) 
which is used in the calibration. A correction for this ef- 
fect would increase the Gamow- Teller strengths extracted 
from ^^Ni(t,^He) (and ■''^Ni(d,^He), since the same pro- 
cedure for the strength calibration was used) by about 
25%. 

Finally, the differences between the various experimen- 
tal and theoretical strength distributions were investi- 
gated in terms of electron-capture rates in the stellar 
environment. At low densities and corresponding nar- 
row electron-capture window, the rates are very sensi- 
tive; to the details of the strength distribution at the 
lowest excitation energies and, therefore, the rates cal- 
culated with the strength distributions predicted us- 
ing the KB3G interaction arc closest to those predicted 
using experimentally determined strength distributions 
from 5®Ni(t,3He) and 5^Ni(d,2He). At higher densities, 
the electron-capture window encompasses a large por- 
tion of the Gamow- Teller strength distribution and the 
rates depend much more on the mean location and the 
width of the Gamow- Teller strength distribution and 
rates calculated using the strength distributions from 
shell- model calculations with the GXPFl or KB3G in- 
teractions do about equally well in reproducing the rates 
calculated with experimental strength distributions from 
the 58Ni(t,3He) and 58Ni(d,2He) data. 
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